Quantum dots in graphene 
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We suggest a way of confining quasiparticles by an external potential in a small region of a 
graphene strip. Transversal electron motion plays a crucial role in this confinement. Properties of 
thus obtained graphene quantum dots are investigated theoretically for different types of the bound- 
ary conditions at the edges of the strip. The (quasi)bound states exist in all systems considered. At 
the same time, the dependence of the conductance on the gate voltage carries an information about 
the shape of the edges. 
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Recently, single layers of carbon atoms (graphene) have 
been obtained experimentally [l|. This new conducting 
material with a high mobility @, HJ has attracted a lot 
of theoretical attention because of its special band struc- 
ture [J, [H, H, 0, H, Q- The spectrum of excitations in 
graphene consists of two conical bands and is described 
by a two dimensional analog of the relativistic Dirac 
equation. 

At the same time, graphene has excellent mechanical 
characteristics and is able to sustain huge electric cur- 
rents [l| . One can attach contacts to submicron graphene 
samples and cut out samples of a desired form and size 
0, 0| . Applying the electric field one can vary consider- 
ably the electron concentration and have both the elec- 
trons and holes as charge carriers. Due to these proper- 
ties the graphene systems look promising for applications 
in nanoelectronic devices. 

One of the most important directions of research us- 
ing the semiconductor heterostructures is fabrication and 
manipulations with so called quantum dots that are con- 
sidered as possible building blocks for a solid state quan- 
tum computer [ioj . The selection of e.g. GaAs/AlGaAs 
for this purpose is related to a possibility of producing 
electrostatic barriers using weak electric fields. Chang- 
ing the field configuration one can change the form of the 
quantum dot, its size and other characteristics. Consid- 
ering applications of the graphene systems, the fabrica- 
tion of the quantum dots looks one the most desirable 
developments in the field. 

In this paper we discuss a method of making quantum 
dots in graphene strips using electrostatic gates. At first 
glance, the possibility of an electrostatic electron confine- 
ment looks surprising since the total density of the con- 
duction electrons is huge n e w 4 x 10 15 cm~ 2 . However, 
the striking feature of the graphene spectrum, namely, 
the existence of the degeneracy points makes the local 
density of the carriers very sensitive to the electric fields. 
This opens a way to create localized states near the zero 
energy of the 2-dimensional Dirac Hamiltonian. 

The existence of bound states in a quantum well is 
one of the basic features of systems described by the 
Schrodinger equation. The situation is different for the 
Dirac equation, since chiral relativistic particles may pen- 



etrate through any high and wide potential barriers. This 
ideal penetration [y, 0] means that one cannot auto- 
matically transfer to graphene the experience in fabrica- 
tion of quantum dots in GaAs using the confinement by 
barriers. 

Fortunately, one can still localize the charge carriers in 
the graphene strip using transversal degrees of freedom 
of their motion. Moreover, in most of the examples below 
the mode ideally propagating along the strip is prohibited 
for the strip of a finite width. 

Formation of the quantum dot in a semiconductor 
wire requires two tunnelling barriers. Surprisingly, in 
graphene it is sufficient to make a single barrier, which 
may be even simpler from the experimental point of view. 
The quasi-bound states exist inside the potential barrier, 
whose left and right slopes work as the " tunnelling barri- 
ers" for the relativistic electrons. The width of the energy 
levels of these quasi-bound states falls off exponentially 
with the width of the barrier and can be very small. 

The very existence of the quasi-bound states (reso- 
nances) is independent of the way of the scattering of the 
electron waves on the edge of the graphene strip (bound- 
ary conditions). However, the positions and the widths 
of the individual resonances and especially the value of 
the background conductance between the resonances de- 
pends on the type of the boundary. Therefore, an exper- 
imental realization of our setup would allow one to study 
properties of the boundary of the real graphene strips. 

The electron wave functions in graphene are usually 
described by a two component (iso)-spinor iji. Its up- 
and down- components correspond to the quantum me- 
chanical amplitudes of finding the particle on one of the 
two sublattices of the hexagonal lattice. In the absence of 
a magnetic field, the usual electron spin does not appear 
in the Hamiltonian and all the electron states have the 
extra double degeneracy. The Fermi level of a neutral 
graphene is pinned near two corners K, K' of the hexag- 
onal Brillouin zone, which generates two valleys in the 
quasiparticle spectrum. The iso-spinor wave function de- 
scribing the low energy electron excitations decomposes 
into a superposition of two waves oscillating with a very 

different wave- vectors tp = e lKr <l>K + e- lK ' r 4>K' 



where 
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FIG. 1: Upper curve: Conductance of the graphene quantum 
dot as a function of Fermi energy for the metallic armchair 
edges for L = 4£, Eq. Q. Lower curve: contribution to con- 
ductance from the transmission channels with p y — ±7r7i/L. 
All calculations are carried out for zero temperature, T = 0. 

4>k = (uk,vk),4>k' = {uk'i Vk') are two smooth en- 
veloping functions. The latter can be found from the 
two-dimensional Dirac equation, e.g. 

[c(p x a x +PyVy) + V(x)](f>K = S(j) K . (1) 

Here c ~ 10 s cm/s is the Fermi velocity and p = —iKV. 
We consider the graphene strip of the width L placed 
along x axis, < y < L. The smooth potential V(x) is 
assumed to be created by an external small size gate (tip). 
We consider the simplest case of a parabolic potential 

V = - {x/xof U/2. (2) 

Details of the asymptotics of the potential at \x\ ^> Xq are 
not important for our results. The envelope function (f>K' 
for the quasiparticle states from the second valley satisfies 
the same Eq. (fT]) with replaced sublattice indices, i.e. 

With <Jy — > — Gy. 

Solution of a couple of two-component Dirac equa- 
tions ([1]) in a strip requires a specification of two bound- 
ary conditions at each edge of the strip. We first consider 
the " armchair" edge corresponding to the boundary con- 
ditions [1| (2/1,2 = 0, L) 

u K \ Vj = e i2 ™iu K ,\ Vj , v K \ yj = e aw 'v K >\ Vi , (3) 

where i = 1, 2 and v\ — 0. If the graphene strip contains 
a multiple of three rows of the hexagons one obtains vi = 
0, which corresponds to a metal. Other numbers of rows 
lead to a semiconducting state with v 2 = ±2/3. 

Eqs. (jl!2p suggest natural units of length and energy 

i= [hcx 2 /U] 1/3 , e Q = he/ (4) 

In the experiments P, d, Q , a graphene strip of the width 
L ps 1/im was separated by a .3/im thick SiC>2 coat- 
ing layer from a n + doped Si wafer. We expect that 
the length scale for the potential V, Eq. ([2]), produced, 
e.g., by varying the thickness of the insulator layer, or 
by local chemical doping, is also xq w lfim. Assuming 
that the characteristic length, Eq. Q, is also £ w 1/im 



we estimate £o w U w .66 x 10~ 3 eV. Making the co- 
ordinate dependent potential, Eq. @, of this strength 
looks rather realistic. For example, reaching the car- 
rier density n s = p^/irh 2 w 10 12 cm~ 2 would require a 
shift of the Fermi energy away from the half-filling by 
AEp — epp ~ .12eV. Even larger carrier densities in 
graphene were reported in the experiments 0, Q ■ 

Solutions of Eqs. (|H2p for the armchair graphene strip 
have a form 

U K - e W»(/ + g) , UK , = e -iPvy/K(f + ( 5 ) 

v K = e ^y' n {f - g) , vk' = e-^y' n {f - g). 

The transverse momentum p y takes the values 

Py( n ) — ( n + u i — v 2 )Ti'hlL , n — ■■, -1, 0, 1, 2, • • ■ . (6) 

The Dirac equation |T]) is now replaced by 

{—ihd/dx + V(x)) f — icp y g — e/, (7) 
(ihd/dx + V{x))g + icp y f = eg. 

These equations decouple from each other and can be 
solved exactly provided the momentum component p y 
perpendicular to the strip vanishes, p y = 0, 

f is -is a f X £ ~ V( x ') 7 / / Q s 
/ = e , g = e , S = / — dx . (8) 

Eq. (j8j) is not what we would like to have because it de- 
scribes the electron waves propagating without reflection 
along the strip. This is just the Id solution considered 
previously [EH, 0]- For p y ^ 0, one cannot solve Eqs. 
(|7|). However, the exact asymptotics at x — ► ±oo of the 
solutions has a simple form {5| 

/- = e lS , <?_ = rer' s , /+ = te? s , g+ = 0, (9) 

where r and t are two complex numbers, |r| 2 + |<| 2 = 1, 
and the subscripts +, — relate to the asymptotics at ±oo. 

The form of the asymptotics chosen in Eq. @ cor- 
responds to the electron flux moving from -oo to +oo, 
where the coefficient r stands for the reflection and t for 
the transmission amplitude. The Landauer formula gives 
the conductance G at zero temperature 

G = G J>„| 2 ,G = 2e 2 //i, (10) 

where t n = t(p y (n)). For the metallic armchair edge, 
the summation goes over n — 0, ±1, ±2, • • • , and |t„| 2 = 
|i_„| 2 . The factor 2 in Go accounts for the electron spin. 

Analytical calculation of the transmission coefficients 
t n is possible only for |e| 3> £o- Fig. 1 shows the depen- 
dence of the conductance on the Fermi energy e calcu- 
lated numerically for L = 4£ for the metallic armchair 
strip. At e > the conductance increases linearly with 
clearly visible steps AG « 2Go corresponding to the 
opening of new channels [TT[ | . At negative Fermi ener- 
gies one can see a series of pronounced resonances. The 
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FIG. 2: Examples of trajectories described in the text drawn 
on the x,p x plane (arbitrary units). Solid lines show the tra- 
jectories with e < either bouncing inside the barrier, or 
reflected by it from the left/right. Tunnelling events between 
the bounded and unbounded trajectories are shown schemat- 
ically (t). Thick dashed lines show the trajectories with e > 
either transmitted for \p y \ < e/c (open channels) or reflected 
for \p y \ > e/c (closed channels). 



resonances appear for all nonzero values of the transverse 
momentum p y ^ © but only those corresponding to 
| n | = 1 and | ii | = 2 are resolved in the figure. 

It is easy to understand the reason for the appearance 
of the (quasi)bound states in graphene. Solutions of the 
equation V[x) — e divide the strip into regions with elec- 
tron or hole type of carriers. The lines separating these 
regions serve as tunnelling barriers for all but normal 
trajectories [f| and this leads to the confinement. 

The (semi)classical dynamics of the massless Dirac 
fcrmions, Eq. ([T]), is given by the effective Hamiltonian 
(see examples of classical trajectories in Fig. 2) 



H cff — e = ±c 



Px + PZ 



V(a 



(11) 



For the (+) sign particles may either fly freely above the 
barrier V(x), or start at the infinity and then be reflected 
from the barrier. The (— ) sign in H c ff corresponds to the 
hole solutions of the Dirac equation, whose trajectories 
bounce inside the barrier for our choice of V(x), Eq. ((2]). 
For a given value of the transverse momentum p y four 
classical turning points where the trajectory changes the 
direction of the propagation along the strip (p x — 0) are 



= ±\ 2 



APv 



= ±\ 2 



-C\p y \ -£ 



U 



(12) 



x V U x 

The electron coming from the infinity is reflected by 
one of the outer turning points a; ou t ± if c \Py\ > e - 
Thus changing the energy e one changes the number 
of open channels, which leads to the conductance quan- 
tization for positive e with the smoothed conductance 
G « 2GoLe/(hnc). 

Finite (hole)trajectories bouncing between the two in- 
ner turning points x- ln± give rise to the quasi-stationary 
states. These trajectories appear at e < — c\p y \ and the 
position of the N-th resonance en may be found from 
the quasiclassical quantization rule 

f " 1+ y/(e N - V{x) f - c*pl- = nh(N + \). (13) 



The resonance acquires a finite width due to quantum 
tunnelling between x; n and io Ut . We can estimate the 
width using a result of Ref. [a], where the transmission 
probability through a linear potential was obtained in 
the form w = cxp(—ircpy/hF), where F = \dU/dx\ is a 
slope of the potential. This result can be used in our case 
provided c\p y \ <§; |ejv|. Since in this case the interval At 
between the reflections at the points x- m _ , x- m+ equals 
At = 2(xq/c) \J— 2£n/U, we find the width 



Tjv = -7— w = 
At 



2xo 



U 



-2a 



N 



■ exp 



TtCP y XQ 

hy/-2e N U 



(14) 



Increasing the characteristic length of the potential Xq 
we can get an extremely narrow levels and long time of 
the confinement of the electrons in such a quantum dot. 

The above results have been obtained for the graphene 
strip with the metallic armchair edges {v% — in Eq. ©). 
Specific for such edges is the existence of the channel with 
p y = providing the perfect transmission at any energy 
\t {e)\ 2 = 1. Below we describe the conductance behavior 
for few other kinds of the edge. 

Fig. 3 shows the conductance of the semiconductor 
armchair graphene strip {y% — ±2/3 in Eq. ([3|)) as com- 
pared to the conductance of the metallic one, both found 
from Eqs. (f7 HT0|> . Several striking differences between 
the two kinds of the edges are clearly seen in the figure. 
First, since there is no channel with p y = 0, the back- 
ground conductance around the resonances at e < van- 
ishes for the semiconductor strip, Gc Go- In the metal- 
lic case the averaged conductance at e < is G w Go- 
Second, the height of the conductance steps at e > is 
AG = 2Go for the metallic graphene strip and AG = Go 
for the semiconductor one. Third, the length of the con- 
ductance plateaus is constant in the metallic case. On 
the contrary, the conductance steps in the semiconduct- 
ing strip have alternating short and long plateaus with 
Ae2 ~ 2Aei. In the experiment, one can expect that the 
metallic and semiconductor strips will be produced in a 
proportion 1:2. 

A way to define the boundary of a Dirac billiard was 
proposed many years ago in Ref. [12j | by introducing 
an infinite mass for quasiparticle behind the boundary. 
Ref. 0] suggested that in graphene this boundary would 
correspond to the transverse confinement of carriers by 
lattice straining. The two (K, K') valleys in this case are 
decoupled from each other and one has [l5[ 

MO) = v K (0) , u K (L) = -v K {L), (15) 



u K + v K = f(x) cos 



PyV 



Uk — Vk = g{x) sin 



PyV 



where the functions / and g are the solutions of Eq. |(7J). 
The boundary conditions (fT5|) are satisfied for p y = 
(n + l/2)nh/L, n = 0, 1, 2, • • • . Each solution, Eq. fl5]). 
is fourfold degenerate. Since p y ^ 0, the conductance 
around the resonances is zero, G <C Gq. The curve G(e) 
is now very similar to what we have found for the metal- 
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FIG. 3: Conductance of the graphene quantum dot for the 
semiconductor armchair edge (solid line) compared to the 
metallic armchair edge (dashed line) for L — 2£ (Q. The 
non-resonant conductance at e < is G « Go for the metallic 
strip and G ~ for the semiconducting one. At e > the 
conductance steps in the semiconductor case are two times 
smaller, AG ~ Go . Alternating series of short and long steps 
are clearly seen for the semiconductor strip. 

lie armchair edges (Fig. 1) provided the picture is shifted 
vertically by —Go- 

Another widely considered type of the edge in graphene 
is the zigzag edge. Since in the case of the zigzag bound- 
ary the edges of the strip belong to the different sublat- 
tices, the components u and v of the envelope function 
vanish at the opposite sides of the strip 

ujc(0),ujp(0) = , v K {L),v K ,{L) = 0. (16) 

In addition to the solutions described by the Dirac equa- 



tion |T]), the zig zag edge supports a band of zero energy 
edge states [13, [H| ■ The (unknown) conductance ~ G 
due to the edge states should be added to the bulk con- 
ductance (|10[) . Except for this edge states contribution 
we do not expect significant differences between the con- 
ductance of the strips with the zigzag edges and those 
confined by the lattice straining [16j. 

In Figs. 1,3 the heights of the resonances are de- 
termined by the level of degeneracy (two- or fourfold) 
of the bound states of the non-interacting electrons. 
In the experiment the shape of the resonances will be 
governed by the electron interaction via the Coulomb 
blockade effect [lj], EH- Since the widths of the reso- 
nances corresponding to large values of the transverse 
momentum (|14[) become exponentially small, the multi- 
ple charging and repopulation of a broad level introduced 
in Ref. [191 ] may occur here. 

To conclude, we considered a possibility of localizing 
charge carriers in a graphene strip by applying an ex- 
ternal electrostatic potential. Such a quantum dot can 
be fabricated using a parabolic potential with a single 
maximum(minimum). Depending on the position of the 
Fermi energy, such a device can serve as either a quan- 
tum dot, or a quantum point contact. The two regimes 
correspond to either the resonance conductance or the 
quantized (step-like) one. An experimental realization of 
our findings would open a way to investigate in graphene 
the reach physics of individually prepared quantum dots. 

This work was supported by the SFB TR 12. Discus- 
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